Nonequilibrium dynamics of a two-channel Kondo system due to a quantum quench 
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Recent experiments by Potok et al. have demonstrated a remarkable tunability between a single- 
channel Fermi liquid fixed point and a two-channel non- Fermi liquid fixed point. Motivated by this 
we study the nonequilibrium dynamics due to a sudden quench of the parameters of a Hamiltonian 
from a single-channel to a two-channel anisotropic Kondo system. We find a distinct difference 
between the long time behavior of local quantities related to the impurity spin as compared to that 
of bulk quantities related to the total (conduction electrons -I- impurity) spin of the system. In 
particular, the local impurity spin and the local spin susceptibility are found to equilibrate, but in a 
very slow power-law fashion which is peculiar to the non- Fermi liquid properties of the Hamiltonian. 
In contrast, we find a lack of equilibration in the two particle expectation values related to the total 
spin of the system. 

PACS numbers: 71.27.-|-a, 72.10.Fk, 5.70.Ln 



The behavior of a local spin coupled to one or more 
independent channels of conduction electrons is a classic 
problem in condensed matter physics. It is well known 
that when the local spin has a size S — 1/2 and is cou- 
pled to only a single channel of conduction electrons, the 
spin is completely screened and the many particle sys- 
tem behaves as a Fermi liquid [1]. In contrast, when the 
local spin is coupled to two or more screening channels, 
one has dramatically different behavior where the spin is 
overscreened, and the system exhibits non-Fermi liquid 
properties [2| . Due to the success in realizing nanostruc- 
tures consisting of a spin coupled to one or more reser- 
voirs, there has been a resurgence of interest in these clas- 
sic systems. The primary focus now is on understanding 
their nonequilibrium properties, such as the effect of cur- 
rent flow [3] and their nonequilibrium time evolution 4] . 

In this paper we will study nonequilibrium dynamics of 
a two-channel Kondo system. We are motivated by recent 
experiments by Potok et al. [3] where by tuning external 
gate voltages a local spin could be effectively coupled to 
a single screening channel or to two independent screen- 
ing channels. Thus as a function of gate voltage, both 
single-channel Fermi liquid physics as well as two chan- 
nel non-Fermi liquid physics was demonstrated on the 
same device. We will study what happens when this ex- 
ternal gate voltage controlling the fine tuning is changed 
rapidly in time from an initial value corresponding to a 
single channel Kondo (ICK) system to a final value cor- 
responding to a two channel Kondo (2CK) system, thus 
inducing nonequilibrium dynamics. The time evolution 
of both single particle expectation values, as well as two 
particle expectation values that exhibit non-Fermi liquid 
behavior in equilibrium will be studied. 

Our system consists of two chiral non-interacting 
fermions that constitute the two channels that inter- 
act with the local spin S = 1/2. We will employ the 
Emery- Kivelson mapping onto an interacting resonant 
level model @ . In the past any dynamics using this map- 



ping, besides addressing other physical situations, was 
studied only at the non-interacting Toulouse point Q. 
In this paper, in order to capture any nontrivial dynam- 
ics of the total spin of the system, we will have to move 
away from the Toulouse point. The Hamiltonian H is: 
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Above, i labels channel while a, /? labels spin index, r, a 
are Pauli matrices and ^ipl^(Ta0'4'ii3 is the spin density op- 
erator for the electrons in the i-th channel, while S = r/2 
is the impurity spin operator. The coupling to the leads 
{t) are time-dependent (in an experiment these may be 
tuned by external gate voltages) Q. When /ii = /12, a 
uniform magnetic field couples equally to both the im- 
purity spin as well as the spins in the leads. When 
/12 = 0,hi 0, the magnetic field couples only to the 
impurity spin. We will assume anisotropic couplings 
Jf(t) = Jf{t) = J^(i) ^ Jf{t). It is convenient to 

define, J, 1 — 



J; 



where SJz,± = at the 2CK fixed point. 

We briefly review the steps involved in map- 
ping the above model onto an interacting reso- 
nant level model Q. One defines the canonically 
conjugate variables [(j)ia{x),Ilji3{y)] — iSijSapS{x — 

written as 



y) 



in terms of 
^a{x) = exp( 



which the 



fermions are 
where 



./V2^ 



This ensures that the 



1^-00 dx'W-iaixl) -I- <l)ia{x) 

same species of fermions anti-commute with each other. 
r\ia are the Klein factors that are necessary to ensure 
anti-commutation between different species of fermions. 
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We choose rji^ = exp{iefj where, 6*^ = O;0f^ = 

^^u; = 7^ (^u + Nit) ; = TT (a^u + + ^2;), 

A'^^iQ being the total number of ia fermions. Defin- 
ing Xia — ^ia — , One changes variables to 2xc — 
Xit + XU + X2t + X21, 2xs = Xit - Xn + X21- - X21, 
2x/ = Xit + XU - X2t - X2;, 2xsf = XiT - XU " 
X2t + X2^- Next one performs the unitary transformation 
H WHU, where U = exp [-^5^x^(0)], followed by a 
refermionization of the Hamiltonian into the fermionic 
fields rft = -iS+; d = ^5" (so that dU - ^ = 5^) 

and ^r,^c,sj,sf{x) = e'^'^'^'^e-^^-^'''^ 1^/2^. In what fol- 
lows we will assume that Jf 2j ^1-L ^re time independent 
and (SJ^ — 0. We set Ji_l = and the only time- 
dependence will be in J2±{t). With this we obtain, 



+h2 I dx^pl{x)^s{x) + {hi - /12) (dU 
-2 (P - ^i;;^) (dU~ i") : V4(0)^s(0) : 
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We will assume that the time dependence of J2±(i) is 
that of a quench, J2±{t) = J±9{t). Thus for t < 
we have a ICK system that is described by an inter- 
acting resonant level model. Whereas for < > 0, the 
Hamiltonian is that of a 2CK system where the cou- 
pling of the resonant level to the reservoir of ipsf fermions 
is via [(d^ - d)V's/(0) + h.c] . Thus in the 2CK 

model effectively only half of the resonant level corre- 
sponding to the Majorana fermion a = —i{d^ — d)/\/2 
couples to the conduction electrons, while the other half 
b = (d^ + d)/ -\/2, does not couple. As was pointed out 
in all non-Fermi liquid behavior stems from this pe- 
culiarity of the resonant level, and as we shall see is also 
responsible for interesting behavior in the dynamics. 

To see this note that immediately after the quench we 
have a highly nonequilibrium system, where any local 
degrees of freedom can relax to the ground state only 
via their coupling to the reservoirs. In the 2CK model 
since only half the local degrees of freedom are coupled, 
local quantities relax very slowly, as we shall show in a 
power law manner. Moreover we find that two particle 
expectation values related to the total (bulk -I- local) spin 
of the system do not relax to their equilibrium values. 

Time evolution of local quantities: We will first 
consider the case when the external magnetic field cou- 
ples only to the local spin, so that hi — h and /i2 — 0. 
We will study the time evolution of the local magnetiza- 
tion, and the local spin susceptibility, the latter in the 



limit /i — > 0. To capture the non- Fermi liquid behav- 
ior of the local susceptibility in an equilibrium 2CK sys- 
tem, it suffices to be at the noninter acting Toulouse point 
= TTVp . Therefore the nonequilibrium dynamics of the 
local quantities will also be studied at the Toulouse point. 
Later while studying the dynamics of the total spin, we 
will have to move away from the Toulouse point so as to 
capture non- Fermi liquid physics ■ We define the 

following Green's functions for the local fermion (spin). 
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Denoting the individual elements of the above matrices as 
Gd^Af Gdi dj 



G 



idt — h 



qR,k obey the equation of motion 
1 
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(9) 
(10) 



where o denotes convolution in time, G^{t,t') = 
[G^{t',t)\ , and E^'-^ are the self-energies due to cou- 

pling to the leads. Defining Fj^ ~ — 



a time- 



dependence of the form J2A_{t) = Ji_0{t) imphes, 
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i + 9{t)e{t') -i9{t') + e{t))\ 
- {e{t') + 9{t)) i + e{t)9{t') ) 



where T is the temperature of the conduction electrons. 

The solutions to Eq. ^ depend on whether the time 
arguments in G^{t,t') are before or after the quench. 
When both times are before the quench, 



G^{t <0,t' <Q)^ 
-^9{t-t')e''-^i'-'') 



(13) 



ih(t—t'^ Q 

e"'(*-*') 



When both times are after the quench we get, 

G^it > 0, i' > 0) = ~i9it - i')e^^(*^*') 
Aiit,t') Bi{t,t') 
Bi{t,t') [Ai{t,t')]* 



(14) 



where Ai{t,t') 



cosh 



ih 



■■ sinh 



Bi{t,t') 
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r^/2 



: sinh 



h^{t-t') 



When one of 



the times is after the quench and the other before, 



G^{t > 0, t' < 0) = -iOit ~ t')e-^+^ 

A2{t,t') B2{t,t') 



(15) 



where [1, [13] x'sioc 2CK — lfT\ 1^ TT equihbrium 
(non-Fermi hquid) local susceptibility of the 2CK Kondo 
system. Thus we find that the logarithmic singularity as- 
sociated with the 2CK system is cut-off by max (^T, ^-^ . 

Moreover at long times TT^ ^ 1, the local susceptibility 
equilibrates, but in a very slow power law fashion which 
is determined by the temperature of the leads. 



cosh 



■■ sinh 



where A2{t,t') 

ih{t-t') 

^B2{t.n= !,f' Binh 

We first discuss the behavior of the magnetization 
5^ = "T'^d'dt(^'0 fo'^ ^ time after the quench (hence. 



t > Q) and when the conduction electrons are at a tem- 
perature T = 0. Substituting Eq. (HH), (HH), (HH) in 
Eq. (Uni), in the fimit of ft, < Tj^ and long times t > 
we find, 



Time evolution of bulk + local quantities : 

Let us consider the case where an external magnetic 
field couples to the total (conduction electrons -f lo- 
cal) spin of the system so that /ii = ft,2 = h. We 

will discuss the time evolution of the response function 

J (LaIY — K^lt — t'Y °^ total spin of the system when ft, — 0. From 

Eq. [31 this may be formally defined as x^i^it'iy^') = 

^idit - t'){{^lix,t)Mx,t).^4iy,t')My,t')})- At the 

Toulouse point = ttvf, the local degrees of freedom 
do not couple to the bulk field V's, so that the response 
function is independent of the local quench and is given 
by the Lindhard function, 



2CK,eq 
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(16) 
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yl-nvp ) qvp — + iS) 



qvp 



(19) 



where 5*^ 



2CK,eq 



In 



rl-2h^+r^y/rl 



-2h'^-T^^jT\-ih^ 



2iT^T\-ih 

is the local magnetization in the ground state of the 2CK 
Hamiltonian. Thus the local magnetization does equili- 
brate, but at a slow rate of h? /T±^ associated with the 
h fermion. In contrast, for a reverse quench J2±(t) = 
J±6{—t) where the time evolution is governed by a ICK 
model for which the a and b fermions are equally coupled 
to the reservoirs, we have checked that Sz equilibrates at 
the much faster rate of rx/4. 

We will now study the time evolution of the local lon- 
gitudinal spin response function xfod^T^') ~ ~~iO{t — 
t'){{S{t)d{t),dHt')d{t')}) which we rewrite as, 

-G^t ,,t it. t')GUt', t) - G% it, t')Gijt', t)] (17) 

It is useful to define the nonequilibrium static susceptibil- 
ity at time T„, xs.iodTm) = d'^xHci^rn + §, " §)• 
For ft = 0, and for very low temperatures T <C T± of the 
conduction electrons, we find the following behavior for 
the static susceptibility at times T,n ^ 1/rj., 



XS,lociTm) Xs.loc,2CK 



-o 



1 



1/1 



(18) 



VTT™ < 1 



ttFi \2TT, 



yTT„, > 1 



Thus the static spin susceptibility xsoiq,^ = 0) = 
' and is independent of q. 



To obtain non-Fermi liquid behavior one has to move 
away from the Toulouse point [3, [l3|, which cou- 
ples ips to the local field, and also introduces non- 
equilibrium dynamics in . Defining, x^(g;i,t') = 
J dx J dy cos g(x — y)x^ix, t; yt')), the leading correction 
in ( — TTVp) to x^ (shown in Fig. [T]) is. 



Ximp(9i t,t ) 



(J- 



t:vf j 



dti / dt2dxdy Xk,citi,t2) 



cos iqix - y)) [G^,^^ ix, 0; t, ti)G^,^^ (0, x; h , t) 
+G'^,^^ix,0;t,h)G^,^^iO,x;h,t) 
G«t^ iO,y;t2,t')G^, iyA.t'M) 



+G^,^^iO,y;t2,t')G^,^^iy,0;t',t2) 



(20) 



where we have assumed that the interaction (J^ — nvp) 
has been switched on adiabatically slowly at long times in 
the past. The label Ximp signifies that it is the correction 
to the bulk response-function due to coupling to the local 
impurity, G^f^ are the Green's functions of the free tps 

fermions, and xfoc defined in Eq. 1171 Defining t = Tm -\- 
^;t' = Tm — § , the nonequilibrium static susceptibility 
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x,t 




0,ti ^ 0,t2 




FIG. 1: Diagram for x?mp{xt,yt'). Dark line: propagator for 
the conduction electrons ^s- Dashed line: propagator for the 
local d fermion. 
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XS,imp{q,Tm) = - [J" -TTVf 

dh I dt2xfoc{tiM) 



(21) 



i^g-2iT„e+ie(ti+t2) 



27r 



where is given in Eq. (jl9p . In equilibrium i.e., when 
J2 = J± and is independent of time, x£c i^ indepen- 
dent of <i + 12. Thus the time integral over ti +t2 forces 
e = in Eq. ([2T|) . With this one recovers the equilib- 

rium result XS^^mM = \ [^^2^)\%\.,2CK = 

Xsim-p 2CK- '^'^ study the evolution of the static suscep- 
tibility after the quench, it is convenient to change vari- 
ables in Eq. ^ io T ^ = <i - ^2- Defining 

u — ^ ) ^^'^ performing the integration over e we 

get, 

,.2 



,R 



Xs 



U 

T 

T 

2 



dT' / dr 



XioAT' + 7:,T'~-) [qvF sin {2qvp\T' - r,„|)](22) 



We will present results for qvp <C T± and times Tm >• 
1 /r_L so that terms that fall off as p-j-^ or faster will be 



dropped. Further, we will consider two cases: one where 
q — 0, and the other when qvp ^ (T, ^r-). 

For q = 0, note that we should first perform the T' 
integral in Eq. (f22|) . and then set q = 0. This gives. 



XS,irnp{q — 0) — cyXs,imp,2CK 



1 



eq 



(23) 



where Xs 



eq 



TnK- — — ■^^f \ _i is the static sus- 
ceptibility in the ICK ground state. For the case qvp ^ 
7^) ' dropping terms of O (y2^j:;^:T~) ' ^"^^ 



Xs,tmpiqvF) = x's. 



1 -i eq 

~^ 2''^ S, imp, ICK 




VII 

q vj. 


SttF^ 


[ n 




q' 
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■In- 



2T 



In^) 
2T J 



2 2\ 9{TTra) 



-Cosint (iTTrn) 
cos {2qVFTrn) 

(24) 



where g{x <^ 1) 1 + 0(a;^), ^(a; > 1) i, and 
. . . represent terms that are small in comparison to 
ln(^),ln(£^). 

Thus we find a marked difference between the suscep- 
tibility at long times after the quench and the suscep- 
tibility in equilibrium Xs%p,2CK- While Xs%p,2CK is 
independent of wave-vector, the out of equilibrium sus- 
ceptibility is strongly dependent on q, and does not even 
reach a time independent steady state, but instead os- 
cillates at frequency qvp (Eq. ([24])). For intermediate 
times TTm ^ 1, performing a time-averaging so that 
terms that oscillate at qvp go to zero, we find. 



XS,imp{qVF > 7^ > T) 



A-kT I 



In- 



1 



AV^T„ 



2 2 

9 «F 



5(ri + g2„^) 



(25) 



Thus for an intermediate time which is longer, the 
lower the temperature, the logarithmic divergences asso- 
ciated with the bulk susceptibility XS,impiqvF) not only 
get cutoff by inverse-time (a result similar to Eq. [18] 
for the local susceptibility), it also acquires some q- 
dependent corrections. In contrast, at long times TT„i ^ 
1, Eq. ([M)) implies that the time-averaged susceptibility 
at large wave-vectors qvp ^ T is, XS,imp{qvF ^ T ^ 

T~) ~ X^simp 2CK^^ (t7^) ' ^"^^ therefore equilibrates. 

The g = static susceptibility (Eq. (^5)) ) on the other 
hand is found to reach a time independent steady state 
which is an equal mixture of the non-analytic in tem- 
perature form of the 2CK ground state, and the ana- 
lytic in temperature form of the ICK ground state. This 
lack of equilibration in bulk properties is consistent with 
nonequilibrium time evolution in integrable models where 
the system retains memory of its initial state. For local 
quantities on the other hand (Eq. ((T6| . dTSH . at least at 
the Toulouse point, the rest of the system to which they 
are coupled acts as a reservoir causing them to equili- 
brate, but at very slow rates compared to a ICK model. 

In summary we have studied the nonequilibrium dy- 
namics in a 2CK system due to a quantum quench. Our 
results highlight how the non-Fermi liquid properties of 
the system, along with its integrability affect the time 
evolution of single particle and two-particle expectation 
values. An interesting question concerns the observabil- 
ity of the nonequilibrium dynamics presented here. Ex- 
periments may be characterized by two kinds of effects 
that have not been taken into account in the present 
treatment. One is that the system could be "open" i.e, 
coupled to some other modes such as phonons, leading 
to an external dissipation rate ^diss which will eventually 
cause the system to equilibrate. The second effect could 
be deviations from integrability arising for example due 
to a nonlinear dispersion for the conduction electrons. 
Studying the consequence of these effects is very inter- 
esting and beyond the scope of this paper. However, one 
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may still be able to speculate on the effect of an exter- 
nal dissipation. In particular a characteristic of the 2CK 
system is slow power-law dynamics. Thus we expect that 
for weak dissipation ^diss ^ T-Li the system will equili- 
brate slowly as / ^ , (where Tm is the time after the 
quench) so that a nonequilibrium/transient state can still 
exist for long enough time-scales to be observable. The 
results of this paper are also relevant for Kondo systems 
in cold-atom gases where dissipative effects are weak [ij). 
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